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ABSTRACT.  A  blast  wave  overpressure-time  history  contains 
voluminous  data  which  can  be  greatly  simplified  by  converting 
them  into  equivalent  yield.  This  calculation  requires  two  in¬ 
dependent  parameters,  usually  chosen  from  measurements  on  peak 
overpressure,  impulse,  standoff  distance,  or  travel  time.  Re¬ 
presentative  computations  are  described.  Concurrence  among  the 
results  of  these  for  a  particular  explosion  indicates  a  simple 
explosion  situation  and  reliable  data;  divergence  implies  some 
effect  such  as  focussed  blast,  or  perhaps  errors  in  measurement. 
All  such  information  is  useful  in  evaluation  of  blast  wave  data, 
and  can  contribute  to  our  understanding  of  explosion  phenomena. 
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EXPLOSIVE  YIELD 


The  magnitude  of  an  explosion  is  described  by  explosive  yield.  A 
representative  example  is  description  of  the  first  nuclear  explosion  as 
equivalent  to  twenty  thousand  tons  of  TNT.  Explosive  yield  serves  to 
specify,  indirectly,  characteristics  of  the  blast  wave  generated  at 
different  distances  from  the  explosion.  Determination  of  explosive 
yield  is  important  in  many  instances;  examples  include  evaluation  of 
some  new  explosive,  assessment  of  hazards  associated  with  manufacture 
or  transportation  of  explosives,  or  development  of  some  new  explosion 
weapon . 


BLAST  WAVE  YIELD 


An  explosive  yield  computed  from  blast  wave  measurements  is  iden¬ 
tified  as  a  blast  wave  yield.  Computation  of  a  blast  wave  yield  is  the 
inverse  of  a  conventional  problem  in  the  physics  of  explosions,  one 
where  blast  wave  characteristics  are  calculated  for  a  specified  explo¬ 
sion.  This  conventional  problem  has  been  solved  theoretically  for  the 
special  case  of  a  centrally  initiated  bare  spherical  charge  exploding  in 
a  uniform  atmosphere,  and  for  related  situations  such  as  a  hemispherical 
charge  exploding  on  an  extended  unyielding  surface  (Ref.  1  and  2). 

These  calculated  values  have  been  adequately  verified  by  experimental 
measurements  (Ref.  3).  These  values  describe  a  reference  explosion,  and 
form  the  basis  for  calculations  of  this  present  report. 

To  convert  values  for  the  reference  explosion  so  that  they  describe 
the  blast  from  some  actual  explosion,  two  independent  parameters  need  be 
specified.  Conventionally  these  are  an  explosive  yield,  usually  ex¬ 
pressed  as  a  multiple  of  the  yield  of  the  reference  explosion,  and  a 
distance  from  the  explosion  center.  Also  needed  are  a  formulation  of 
the  scaling  law  for  explosions ,  plus  a  complete  mathematical  description 
of  the  reference  explosion.  Together  these  items  establish  the  entire 
pressure-time  history  for  the  particular  blast  wave  system  of  concern. 
Details  of  the  computation  are  described  in  several  sources,  among  them 
Ref.  3,  4,  and  5. 
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THREE  TYPES  OF  BLAST  OVERPRESSURE 


In  description  of  a  particular  blast  wave  system,  the  blast  pres¬ 
sure  item  may  be  described  in  terms  of  either  an  absolute  pressure  or  as 
an  overpressure,  the  pressure  above  ambient.  In  this  material  we  work 
exclusively  in  terms  of  blast  overpressure,  an  item  which  in  practical 
work  may  be  referred  to  by  the  simple  term  "pressure."  We  should  note, 
however,  that  there  are  three  different  types  of  blast  wave  overpressures 
corresponding  to  three  different  arrangements  of  the  gages  used  for  the 
measurement . 

Simplest  of  these  blast  pressures  is  the  free  field  value  that 
exists  in  the  undisturbed  blast  wave.  A  gage  with  sensing  element 
normal  to  the  direction  of  propagation  of  the  blast  (side-on  to  the 
blast  front)  and  which  also  does  not  disturb  the  blast  system  can  mea¬ 
sure  this.  Thus,  the  phrase  "side-on"  commonly  indicates  a  free  field 
value.  Results  from  a  pencil  type  of  blast  gage  offer  an  example.  Here 
the  physical  arrangement  of  the  gage  together  with  its  small  size  and 
pointed  front  give  measured  side-on  values  that  can  be  very  close  f<-> 
those  of  the  free  field. 

Other  blast  gage  arrangements  involve  face-on  gages.  Here  the 
sensing  element  faces  the  explosion  center  and  the  blast  impinges  di¬ 
rectly  onto  it.  However,  depending  on  the  physical  arrangement,  there 
are  two  different  face-on  overpressures.  One  of  these  has  a  sensing 
element  that  is  part  of  an  extended  unyielding  surface  which  also  is 
normal  to  the  direction  of  propagation  of  the  blast.  The  large  un¬ 
yielding  surface  causes  the  blast  wave  to  be  reflected  so  that  the 
measurement  is  of  the  reflected  overpressure.  A  face-on  gage  at  a  hard 
earthen  surface,  or  perhaps  one  embedded  in  a  large  sheet  of  armor 
plate,  can  provide  a  reflected  overpressure-time  history  for  an  imping¬ 
ing  blast  wave  system. 

Another  type  of  measurement  uses  a  face-on  sensing  element  without 
the  surrounding  reflecting  surface.  Here  the  gage  acts  in  the  manner  of 
one  of  the  elements  of  a  Pitot  tube,  and  measures  the  sum  of  the  free 
field  overpressure  and  an  impinging  effect  of  the  blast  wind.  This  is 
termed  a  stagnation  overpressure,  implying  that  the  measuring  device 
has  rendered  the  blast  wind  stagnate.  A  stagnation  overpressure  is 
necessarily  greater  than  the  free  field  overpressure,  and  is  necessar¬ 
ily  less  than  the  reflected  overpressure. 


BLAST  MEASUREMENT  PROBLEMS 


Computation  of  a  blast  wave  yield  requires  data  such  as  the  over¬ 
pressures  generated  at  a  particular  location.  All  the  various  types  of 
gages  for  this  measurement  have  special  problems  of  their  own,  and  there 
are  associated  problems  in  recording  instruments  and  interconnecting  de¬ 
vices.  One  problem  for  the  gage-connection-recorder  system  is  response 
time,  for  the  system  must  register  pressures  nearly  instantaneously,  yet 
be  without  overshoot.  Calibration  of  the  system  is  also  something  of  a 
problem,  perhaps  more  severe  than  sometimes  realized.  Calibration  of 
the  system  is  necessary  in  order  to  convert  its  output  into  overpres¬ 
sures,  and  ideally  this  calibration  should  be  made  on  a  time  scale  com¬ 
parable  with  that  of  the  desired  measurement.  Yet  it  is  difficult  to 
obtain  pressure  rise  rates  in  a  calibration  device  that  are  at  all  com¬ 
parable  with  those  of  an  explosion. 

Subtle  difficulties  with  blast  measurements  .re  also  associated 
with  details  of  gage  placement.  Thus  gage  orientation  and  type  of  over¬ 
pressure  being  measured,  free  field,  reflected,  or  stagnation,  must  be 
specified  for  the  data  to  be  meaningful.  A  gage  with  a  side-on  sensing 
element  should  not  interact  with  the  blast  but  may  actually  do  so  and 
cause  observed  values  to  differ  from  those  of  the  free  field.  One 
troublesome  source  of  possible  interaction  is  the  mounting  for  the  gage, 
which,  of  course,  must  be  sturdy,  but  not  sufficiently  large  to  cause  a 
disturbance.  Face-on  gages  are  also  subject  to  special  difficulties. 
Thus,  a  gage  may  protrude  from  a  surrounding  surface  so  that  the  type 
of  overpressure  it  measures  changes  from  stagnation  to  reflection  dur¬ 
ing  the  measurement.  Alternatively,  a  reflecting  surface  surrounding 
a  gage  may  be  too  small  to  block  all  flow;  here  the  type  of  overpressure 
being  measured  changes  during  the  measurement  from  that  of  a  reflected 
to  a  stagnation  overpressure.  In  such  cases  interpretation  of  the  d?ta 
is  difficult. 

Blast  gages  also  are  particularly  sensitive  to  the  presence  of 
neighboring  surfaces,  for  these  can  give  unwanted  reflections  of  annoy¬ 
ing  nature.  This  aspect  may  place  limitations  on  the  physical  arrange¬ 
ment  for  a  test,  as,  for  example,  in  the  vertical  height  of  test  gear 
or  the  location  of  support  structures.  Furthermore,  some  blast  gages 
with  associated  connections  are  sensitive  to  acceleration;  here  motion 
caused  by  the  blast  can  introduce  error.  In  view  of  all  these  and  other 
difficulties,  it  is  only  to  be  expected  that  blast  overpressure -time 
histories  ordinarily  are  subject  to  uncertainties.  Some  aspects  of  this 
situation  are  considered  below. 

A  typical  blast  overpressure-time  record  and  its  required  calibra¬ 
tion  are  illustrated  in  Fig.  1. 
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FIG.  1.  Typical  Overpressure-Time  History 
and  Pressure  Calibration. 
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OVERPRESSURE-TIME  HISTORY 


An  overpressure-time  record  for  a  blast  wave  contains  a  great  deal 
of  information  about  the  explosion.  The  problem  here  is  to  organize 
this  information  into  simple  meaningful  terms.  For  this  it  is  conven¬ 
ient  to  extract  from  the  overall  record  two  pertinent  values,  (1)  the 
peak  overpressure  that  exists  at  the  blast  wave  front,  and  (2)  the 
overpressure-time  integral  for  the  initial  positive  phase  of  the  blast 
wave  system.  This  latter  represents  a  positive  impulse  per  unit  area, 
and  may  be  referred  to  by  the  simple  term  "impulse."  Other  items  such 
as  duration  of  the  positive  pressure  phase  and  its  rate  of  decay  may 
also  be  of  concern. 

The  peak  overpressure  for  the  blast  wave  can  be  found  by  back 
extrapolation  of  the  initial  portion  of  the  record  to  time  zero,  the 
instant  of  arrival  of  the  blast  wave  front  (Ref.  6).  This  extrapola¬ 
tion  avoids  difficulty  with  overshoot  or  with  finite  rise  time  in  the 
gage-connector-recorder  system.  For  this  purpose  it  is  desirable  to 
convert  the  analogue  record  as  illustrated  above  to  digital  form,  and 
then  perform  the  extrapolation  by  computer.  Figure  2  illustrates  this 
conversion.  Here  readings  taken  from  the  analogue  record  and  its  cali¬ 
bration  have  been  converted  to  overpressure  and  time,  and  plotted  as 
points  on  semilogarithmic  coordinates  for  the  first  two-thirds  of  the 
record.  A  straight  line  is  then  fitted  by  a  least-squares  technique. 

The  intercept  of  this  line  gives  the  peak  overpressure;  its  slope  de¬ 
scribes  the  decay  characteristic  of  the  blast  wave.  The  intercept  in 
this  case  indicates  a  peak  overpressure  of  58.3  psi,  versus  a  maximum 
observed  overpressure  of  38.5  psi,  and  the  slope,  m,  is  -2.25  ms-1. 

The  discrepancy  between  peak  overpressure  and  maximum  overpressure  as 
illustrated  here  is  typical  of  a  situation  where  a  low  impedance  rela¬ 
tively  slow  gage  system  is  used.  This  discrepancy  indicates  some  of 
the  possible  errors  associated  with  interpretation  of  blast  wave  data. 

A  parallel  technique  using  the  last  third  of  the  positive  over-' 

,  res sure  phase  permits  a  precise  identification  of  blast  wave  duration, 
tj •  This  is  illustrated  in  Fig.  3,  where  a  duration  is  found  as  1.39  ms. 
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FIG.  3.  Duration  Time  Found  by  Smoothing. 
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From  results  of  such  calculations  the  record  for  the  entire  posi¬ 
tive  phase  of  a  classical  overpressure -time  history  can  be  fitted  to 
the  semi -empirical  equation  (Ref.  3) 


P 


p°(l-t/td)e 


-bt/t 


(1) 


where  p  is  the  overpressure  at  time  t  from  zero  time,  pu  the  peak  over¬ 
pressure  that  occurs  at  zero  time,  t.^  the  duration  of  the  positive  phase 
and  e  the  base  of  natural  logarithms.  Here,  the  item  b  is  a  decay  pa¬ 
rameter  related  to  the  slope  m  of  a  previous  figure  and  found  as  b  = 

-  mt^  -  1,  and  equals  2.14  for  the  blast  wave  above.  Figure  4  is  recon¬ 
structed  from  the  data  of  Fig.  1.  This  shows  observed  values  as  points, 
and  a  line  that  corresponds  to  Eq .  1.  General  concurrence  is  evident. 
Validity  of  the  equation  and  of  its  defining  parameters  can  be  further 
checked  by  values  for  impulse,  that  is,  the  overpressure-time  integral. 
For  this  one  compares  the  results  of  an  analytic  integration  of  the 
equation  with  those  obtained  by  numeric  integration  using  the  original 
data.  In  this  instance  the  impulse  indicated  by  integration  of  Eq.  1 
is  22.3  psi-ms,  and  that  found  by  direct  numeric  integration  is  21.4 
psi-ms.  Reasonable  agreement  such  as  this  is  to  be  expected  in  most 
circumstances  where  the  blast  wave  configuration  is  classical  in  nature. 
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THEORETICAL  ASPECTS  OF  BLAST  WAVE  YIELD  COMPUTATION 


The  two  items  obtained  above,  (1)  peak  overpressure,  and  (2)  im¬ 
pulse  are  important  for  many  purposes,  as  for  example,  in  assessment  of 
damage  potential.  In  experimental  studies  on  blast  two  additional 
items,  (3)  distance  from  center  of  explosion  to  sensing  element,  or 
standoff  distance,  and  (A)  the  time  interval  between  initiation  of  the 
explosion  and  arrival  of  the  blast  front  at  the  gage,  identified  as 
travel  time  or  as  time  of  arrival,  frequently  are  readily  measured  and 
provide  additional  useful  information.  An  important  observation  con¬ 
cerning  the  four  items  is  that  only  two  of  them  are  needed  to  establish 
a  value  for  blast  wave  yield.  This  observation  is  consistent  with  item 
(1)  above,  that  two  independent  parameters  can  specify  a  blast  wave.  Now 
there  are  six  possible  combinations  of  four  things  taken  two  at  a  time. 
Hence,  the  four  items  here  can  provide  six  separate  estimates  for  the 
blast  wave  yield.  If  more  than  one  gage  is  used  for  the  test,  then, 
more  than  one  set  of  six  estimates  can  be  made. 

A  useful  aspect  of  such  multiple  estimates  of  blast  wave  yield  is 
that  they  permit  an  immediate  evaluation  of  the  internal  consistency  of 
observed  data.  If  the  estimates  agree  reasonably  well  with  each  other, 
confidence  can  be  placed  in  the  overall  results.  On  the  other  hand, 
disagreement  may  indicate  errors  in  measurements  or  in  calibrations. 

It  also  may  indicate  a  different  type  of  explosive,  or  blast  anomalies, 
such  as  jetting  or  an  asymmetrical  wave  front,  or  that  the  configuration 
of  the  actual  blast  wave  system  differs  from  the  classical  one  of  the 
reference  explosion.  Disagreement  may  result  from  effects  such  as  a 
blast  being  focussed  to  increase  its  peak  overpressure  without  a  cor¬ 
responding  increase  in  impulse,  or  the  presence  of  inert  material  in 
the  explosive  that  increases  reflected  impulse  without  a  corresponding 
increase  in  side-on  overpressure  . 

Another  aspect  of  these  blast  wave  yield  calculations  of  possible 
interest  concerns  distributed  energy  explosions.  Here  there  may  be  no 
definite  explosion  center  and  no  definite  travel  time.  However,  over¬ 
pressure-time  history  above,  through  peak  overpressure  and  impulse, 
makes  it  quite  feasible  to  find  an  effective  yield.  Then  by  inverting 
the  yield  calculation  an  effective  standoff  distance  can  be  determined. 
This  is  illustrated  by  a  calculation  in  the  Appendix.  Thus,  in  many 
circumstances,  determination  of  blast  wave  yields  can  contribute  both 
to  an  appreciation  of  the  instrumentation  problems  in  explosion  tests 
and  to  our  understanding  of  explosions  and  their  damage  potential  in 
general  . 
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SCALING  LAW  EXPRESSIONS 


Determination  of  explosive  yield  from  blast  wave  measurements  is 
essentially  an  application  of  the  scaling  law  for  explosions.  This  law 
as  used  here  expresses  the  observation  that  an  explosion  creates  a  dis¬ 
turbance  by  distributing  energy  throughout  a  surrounding  volume,  and 
that  two  disturbances  are  similar  if  they  have  the  same  geometry  and  if 
the  energy  release  into  unit  amount  of  disturbed  material  is  the  same. 
Now  for  an  explosion  in  a  uniform  atmosphere  the  total  amount  of  dis¬ 
turbed  material  is  proportional  to  its  volume  (cube  of  a  distance)  and 
also  to  its  density.  Thus,  for  two  similar  but  different  explosion 
disturbances,  identified  by  subscripts  1  and  2, 


Wj /W2  =  p ir?/pj>r$  (2) 

where  W  is  the  explosive  energy  release,  r  the  distance  from  center  of 
the  explosion  to  edge  of  the  disturbance,  (this  distance  is  a  function 
of  time),  and  p  is  the  density  of  the  undisturbed  surrounding  atmosphere. 
This  formulation  of  the  scaling  law  is  a  modification  of  an  earlier  form 
(Ref.  7)  that  uses  pressure  rather  than  density  as  an  index  of  mass. 

Now  consider  a  reference  explosion  with  energy  release  into  an 
atmosphere  with  density  pg.  Let  a  distance  from  this  reference  explo¬ 
sion  be  termed  a  "scaled  distance."  Equation  2  now  indicates  that 
another  explosion  with  energy  release  W  creates  a  similar  disturbance 
in  an  atmosphere  with  density  p  under  conditions  such  that 

scaled  distance  =  (actual  distance)  (p /p g) 1 ^ 3 (W0/W) 1 (3) 

This  important  expression  is  the  basis  for  the  blast  wave  yield  calcula¬ 
tions  described  here.  It  may  be  shown  (Ref.  5)  that  it  contains  within 
it  a  corresponding  definition  of  a  scaled  time,  that  is,  a  time  interval 
pertaining  to  the  reference  explosion,  as 

scaled  time  =  (actual  time)  (p/pg)1  3(W0/W)1  3(a/a0)  (A) 


where  a  is  the  speed  of  sound,  and  the  subscript  0  identifies  a  value 
for  the  reference  explosion.  A  corresponding  expression  for  scaled  im¬ 
pulse  becomes 

scaled  impulse  =  (actual  impulse) 

(p/p0)l/a(Wo/W)l/3(a/a0)(P0/P)  (5) 
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where  P  indicates  absolute  atmospheric  pressure. 

These  forms  of  the  scaling  relations  are  the  inverse  of  conventional 
ones  that  describe  an  actual  explosion  in  terms  of  values  computed  for  a 
reference  explosion. 


YIELD  EXPRESSIONS 


Each  of  the  last  three  formulations  of  the  scaling  law  above  can  be 
rearranged  to  provide  a  direct  solution  of  an  actual  energy  release 
(yield)  W  in  terms  of  a  reference  yield  W0 .  For  this  it  is  convenient 
to  express  the  atmospheric  density  ratio  in  terms  of  absolute  pressure 
and  temperature,  or  (p/p0)  =  (P/P0)Tq/T).  Also  the  speed  of  sound  in 
the  atmosphere  varies  as  the  square  root  of  its  absolute  temperature,  or 
(a/ag)  =  (T/To)1/2.  Three  expressions  for  blast  wave  yield  W  then  be¬ 
come 


W  =  Wff  (actual  distance/scaled  distance) 3 (P/PB ) (T0 /T)  (6) 
W  =  W0  (actual  time/scaled  time) 3 (P/Pg) (T/Ty) 1 ^  (7) 
W  =*  W0(actual  impulse/scaled  impulse)3(P0/P)2  (T/Tg)1^2  (8) 


These  equations  pertain  to  geometrically  similar  blast  wave  systems  that 
create  the  same  intensity  of  disturbances  in  their  surrounding  atmo¬ 
spheres  . 

The  pressure  and  temperature  ratios  of  the  above  three  equations 
can  be  combined  with  measured  values  of  distance,  time,  or  impulse  to 
give  terms  which  are  corrected  for  atmospheric  conditions  if  different 
from  those  for  the  reference  atmosphere.  That  is, 

corrected  distance  =  actual  distance (P/PB ) 1  ^ 3 (Tg /T) ^3  (9) 

corrected  time  =  actual  time(P/P0)  1//3(T/Tb)  ^  b  (10) 

corrected  impulse  =  actual  impulse (Pg /P)2 ^ 3 (T/Tg) 1 ^ 6  (11) 
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These  corrected  items  can  make  an  apparent  simplification  in  the  equa¬ 
tions  for  yield  calculation.  Working  equations  for  blast  wave  yield  W 
in  terms  of  a  reference  yield  Wq  so  become 


W  =  W0 (corrected  distance/scaled  distance)3 
W  =  Wq (corrected  time/scaled  time)3 
W  =*  W,j  (corrected  impulse/scaled  impulse) 3 


(12) 

(13) 

(1A) 


As  a  check  on  these  equations,  note  that  working  Eq .  12  when  combined 
with  Eq.  9  gives  the  basic  scaling  law  of  Eq .  3,  etc. 

The  fractional  power  terms  involved  in  the  corrected  items  of  the 
above  equations  may  in  some  cases  be  unimportant  and  then  can  be  neg¬ 
lected.  It  may  also  be  noted  that  each  blast  wave  yield  as  computed  in¬ 
volves  the  cube  of  a  ratio,  hence  any  uncertainty  in  the  original  mea¬ 
surement  is  magnified  by  a  factor  of  three  in  this  computation.  Blast 
wave  yields,  therefore,  are  necessarily  rather  imprecise. 


REFERENCE  EXPLOSION  DATA 


The  above  equations  are  in  terms  of  a  corrected  item  that  involves 
measurements  on  an  actual  explosion  and  a  scaled  item  for  the  geometri¬ 
cally  similar  blast  wave  of  a  reference  explosion.  Thus,  to  use  these 
relations  a  reasonably  complete  description  of  a  reference  explosion  is 
needed.  This  description  may  be  in  the  form  of  graphs,  tables,  or  even 
equations;  many  such  data  are  scattered  through  the  technical  literature. 
A  fairly  detailed  and  convenient  source  of  data  is  provided  by  the  tables 
of  Ref.  3,  and  these  tables  are  used  in  the  calculations  illustrated 
here.  These  are  in  English  units  with  the  explosion  of  a  bare  spherical 
charge  of  one  pound  of  TNT  in  a  uniform  atmosphere  at  13.6  psia  and  59F 
as  the  reference.  A  second  table  is  available  for  use  with  forthcoming 
metric  units;  this  describes  the  explosion  of  one  kilogram  of  TNT  in  an 
atmosphere  at  one  bar  (1,000  mb)  and  16C  (59F) .  The  reference  tables 
may  be  used  directly  for  manual  computation,  as  illustrated  here,  or 
they  may  be  stored  in  equivalent  condensed  form  for  ready  computer  use. 
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YIELD  CALCULATED  FROM  OVERPRESSURE  AND  DISTANCE 


The  ratio  of  peak  overpressure  to  ambient  pressure  is  a  direct 
measure  of  blast  wave  intensity,  hence  the  scaling  law  requirement  for 
similarity  is  met  where  two  different  disturbances  show  the  same  over¬ 
pressure  ratio.  This  requirement  for  identical  ratios  indicates  that  a 
correction  must  be  applied  if  an  actual  ambient  pressure  differs  from 
that  assumed  for  the  reference  explosion.  That  is, 

corrected  overpressure  =  actual  overpressure  (Pg/P)  (15) 


where  Pg  is  the  absolute  atmospheric  pressure  for  the  reference  explo¬ 
sion,  and  P  the  ambient  pressure  in  the  atmosphere  of  the  actual 
explosion . 

A  yield  calculation  based  on  overpressure  requires,  in  addition, 
data  on  some  other  item,  in  this  case  taken  as  a  standoff  distance. 

This,  however,  must  be  the  actual  distance  from  the  center  of  the  explo¬ 
sion  to  the  sensing  element,  and  is  not  quite  identical  with  distance 
from  charge  surface  to  tip  of  a  gage.  A  convenient  routine  for  computa¬ 
tion  of  yield  using  overpressure-distance  data  is  as  follows.  Note  that 
the  routine  requires  incidental  data  on  the  pressure  and  temperature  of 
the  actual  atmosphere. 

1.  Find  the  peak  overpressure  from  the  blast  wave  record, 
as  by  the  extrapolation  method  described  above. 

2.  Convert  this  peak  overpressure  to  reference  conditions 
by  multiplying  by  the  factor  (Pg/P),  as  in  Eq .  15. 

3.  Enter  the  table  for  the  reference  explosion  at  the 
corrected  value  of  the  peak  overpressure  and  find  a 
corresponding  standoff  distance  for  the  reference 
explosion.  This  is  the  scaled  distance. 

4.  Correct  the  actual  distance  for  atmospheric  pressure 
and  temperature,  as  by  Eq .  9.  This  gives  the  corrected 
distance . 

5.  Compute  yield  as  the  cube  of  the  ratio  (corrected  dis¬ 
tance/scaled  distance),  as  by  Eq .  12. 


The  reference  overpressure  used  in  this  calculation  must  be  of  the 
same  type  as  that  for  the  actual  explosion,  that  is,  free  field  re¬ 
flected,  or  stagnation.  Both  free  field  and  reflected  values  are  pro¬ 
vided  directly  in  reference  tables.  For  the  less  frequently  used 
stagnation  overpressure  the  most  convenient  procedure  seems  to  be  to 
convert  to  free  field  values.  Formulas  and  tables  that  provide  for 
this  conversion  are  given  in  Ref.  5. 
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The  five  steps  of  the  calculation  above  are  illustrated  using 
data  for  the  explosion  whose  overpressure-time  history  is  shown  in 
Fig.  1  through  4,  and  for  which  other  data  are  given  in  the  Appendix. 

Step  1.  The  determination  of  peak  overpressure  from  the  over¬ 
pressure-time  history  is  illustrated  in  Fig.  2.  This  peak 
overpressure  is  found  to  be  58.3  psi. 

Step  2.  By  Eq.  9,  correct  peak  overpressure  =  58.3  *  (13.6/13.6) 
=  58.3  psi  . 

Step  3.  Entering  the  table  of  Ref.  3  at  a  peak  reflected 
overpressure  of  58.3  psi,  the  corresponding  scaled  distance 
is  seen  to  be  between  5.60  and  5.65  feet.  Linear  interpola¬ 
tion,  adequate  for  purposes  here,  indicates  a  scaled  dis¬ 
tance  of  5.64  feet. 

Step  4.  Corrected  standoff  distance  =  5.50  x  (13 .6/13 .6) 1 ^3 
(59+460) /(102+460)1/ 3  =5.36  feet. 

Step  5.  Yield  =  (5. 36/5. 64) 3  =  0.86  pounds  TNT. 


YIELD  CALCULATED  FROM  OVERPRESSURE  AND  TRAVEL 
TIME  OR  IMPULSE 


A  yield  calculation  using  overpressure  with  travel  time  rather  than 
with  standoff  distance  is  quite  similar  to  the  calculation  above,  except 
of  course,  that  Eq .  10  and  13  involving  time  are  used  rather  than  those 
involving  distance.  Calculation  of  yield  from  overpressure  and  impulse 
also  proceeds  in  a  similar  manner.  The  Appendix  includes  the  results  of 
these  plus  other  calculations  of  blast  wave  yield. 


YIELD  CALCULATED  FROM  STANDOFF  DISTANCE 
AND  TRAVEL  TIME  ONLY 


Yield  calculation  using  overpressure  are  straightforward  in  that 
they  permit  direct  entry  into  a  reference  table.  Where  only  standoff 
distance  and  travel  time  are  to  be  the  basis  for  a  yield  calculation, 
such  direct  entry  may  not  be  possible.  This  difficulty  can  readily  be 
circumvented,  however.  Thus,  by  combining  working  Eq .  12  and  13  it  can 
be  seen  that  the  ratio  (corrected  distance/corrected  time)  is  identical 
with  the  ratio  (scaled  distance/scaled  time).  This  ratio  has  the  physi¬ 
cal  meaning  of  average  travel  speed,  and  also  is  a  measure  of  explosion 
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disturbance  intensity.  A  technique  for  computation  of  yield  using  such 
data  is  as  follows. 

Step  1.  Correct  standoff  distance  for  ambient  pressure  and 
temperature,  as  by  Eq .  9 ,  and  travel  time  as  by  Eq.  10. 

Divide  the  first  by  the  second  to  find  the  corrected  (dis¬ 
tance/time)  ratio.  This  is  also  the  scaled  ratio. 

Step  2 .  Enter  the  table  describing  the  reference  explosion 
at  the  computed  (dist ance/time)  ratio,  and  find  a  correspond¬ 
ing  scaled  distance. 

Step  3.  Compute  blast  wave  yield  as  the  cube  of  the  ratio 
(corrected  distance/scaled  distance),  as  by  Eq .  12. 

An  alternative  procedure  is  to  find  scaled  time  in  Step  2  rather 
than  scaled  distance,  and  then  compute  yield  as  the  cube  of  a  time  ratio. 
The  two  procedures  give  identical  results. 

A  convenient  circumstance  for  the  above  calculation  is  that  the 
tables  of  Ref.  3  provide  directly  the  needed  values  for  the  average 
travel  speed  for  the  blast  wave  of  the  reference  explosion.  If  tables 
without  this  item  are  to  be  used,  it  can  readily  be  provided  by  simple 
division  of  tabulated  distance  by  tabulated  time. 


YIELD  CALCULATED  FROM  IMPULSE-DISTANCE  AND 
FROM  IMPULSE -TIME  DATA 


Computing  a  yield  using  only  impulse-standoff  distance  data  has  a 
difficulty  mentioned  above;  there  may  be  no  means  for  ready  entry  to 
the  table  that  describes  the  reference  explosion.  However,  by  combin¬ 
ing  Eq.  12  and  14  it  can  be  seen  that  the  ratio  (corrected  impulse/ 
corrected  distance)  is  identical  with  the  ratio  (scaled  impulse/scaled 
distance) .  A  computational  technique  based  on  this  observation  is  as 
follows . 

Step  1 .  Correct  impulse  and  standoff  distance  for  atmo¬ 
spheric  pressure  and  temperature,  as  by  Eq .  9  and  11.  Then 
find  the  ratio  of  these.  This  also  is  the  ratio  (scaled 
impulse/scaled  time). 

Step  2 .  From  individual  tabulated  values  for  impulse  and 
for  distance  in  the  table  for  the  reference  explosion,  pre¬ 
pare  a  supplemental  column  of  (impulse/distance)  ratios. 

Enter  this  column  with  the  ratio  of  corrected  values  found 
in  Step  1,  and  note  the  corresponding  scaled  distance. 
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Step  3.  Compute  yield  as  cube  of  the  ratio  (corrected  dis¬ 
tance/scaled  distance)  . 

A  parallel  routine  using  data  for  impulse  and  for  travel  time  per¬ 
mits  a  separate  yield  computation,  for  by  Eq .  13  and  14  the  ratios 
(corrected  impulse/corrected  travel  time)  and  (scaled  impulse/scaled 
travel  time)  are  identities. 

In  all  these  calculations  it  is  important  to  know  the  type  of  im¬ 
pulse  measured,  whether  free  field,  reflected,  or  stagnation,  and  to  use 
corresponding  values  for  the  reference  explosion.  It  should  also  be 
noted  that  relatively  few  descriptions  of  the  blast  wave  for  a  reference 
explosion  include  all  these  individual  impulse  items. 


COMPUTER  CALCULATIONS 


The  above  yield  computations,  including  extraction  of  peak  over¬ 
pressure  and  determination  of  observed  impulse,  can  all  readily  be 
performed  by  computer.  A  program  using  digitalized  data  in  punched -card 
form  has  been  prepared  for  this  purpose,  and  its  use  is  indicated  for 
routine  or  repetitive  calculations . 
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Appendix 


TYPICAL  CALCULATIONS 


.  MEASURED  CONDITIONS  OF  TEST  EXPLOSION 

ambient  pressure 
ambient  temperature 
standoff  distance 
travel  time 


13.6  psia 
102  F 
5.50  ft 
1 .71  ms 


Obtained  From  Overpressure -Time  History  (Reflected  Values) 

peak  overpressure  58.3  psi 

impulse  21  .4  psi-ms 


Corrected  to  Reference  Conditions  of  59F  and  13.6  psia 


corrected  peak  overpressure 
corrected  standoff  distance 
corrected  travel  time 
corrected  impulse 


58.3  psi 
5.36  ft 
1 .69  ms 
21.1  psi-ms 


Computed  Yields  (pounds  TNT)  based  on: 

peak  overpressure  and  standoff  distance 
peak  overpressure  and  travel  time 
peak  overpressure  and  impulse 
standoff  distance  and  travel  time 
impulse  and  standoff  distance 
impulse  and  travel,  time 

average 

standard  deviation 


0.86 

0.83 

1.07 

0.86 

0.92 

0.98 

0.92 

±0.06 


The  concurrence  illustrated  is  representative  of  reasonably  good 
measurements  made  under  field  conditions. 
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2.  ANALYSIS  OF  BLAST  FROM  A  FUEL-AIR  EXPLOSION 

The  overpressure-timt.  record  made  at  the  surface  of  the  earth 
during  an  overhead  fuel-air  explosion  indicated  a  peak  reflected  over-  . 
pressure  of  245  psi  and  a  positive  reflected  impulse  of  210  psi-ms. 
Assume  atmospheric  pressure  and  temperature  close  to  those  of  the  refer¬ 
ence  explosion  so  that  correction  may  be  neglected  and  find  the  indi¬ 
cated  yield  and  the  equivalent  standoff  distance. 

This  peak  reflected  overpressure  corresponds,  by  the  tables  of  Ref. 

3,  to  a  scaled  reflected  impulse  of  42.3  psi-ms.  Then  by  Eq.  14, 

indicated  yield  =  (210/42. 3) 3  =  130  pounds  TNT 

This  observed  peak  overpressure  also  indicates  a  scaled  distance 
of  3.65  feet.  By  Eq .  3,  when  inverted 

actual  distance  =>  3.65  (lSO)1^3  =*  18.5  feet 

That  is,  this  distributed  energy  explosion  is  equivalent  at  this 
particular  location  to  the  concentrated  energy  explosion  of  a  bare 
spherical  charge  of  130  pounds  of  TNT  directly  overhead  at  a  height  of 
18.5  feet. 


NWC  TP  5020 


REFERENCES 


1.  Brode,  H.  L.  "Numerical  Solution  of  Spherical  Shock  Waves,"  J  APPL 
PHYS,  Vol .  26  (1955),  p.  766. 

2.  Rand  Corporation.  A  Calculation  of  the  Blast  Wave  From  a  Spherical 
Charge  of  TNT,  by  H.  L.  Brode.  Santa  Monica,  Calif.,  RAND,  August 
1957.  (Research  Memorandum  RM-1965.) 

3.  Naval  Weapons  Center.  Engineering  Elements  of  Explosions,  by  G.  F. 
Kinney.  China  Lake,  Calif.,  NWC,  November  1968.  (NWC  TP  4654.) 

4.  Glasstone,  S.  The  Effects  of  Nuclear  Weapons.  U.  S.  Atomic  Energy 
Commission,  1962. 

5.  Kinney,  G.  F.  Explosive'  Shocks  in  Air.  New  York,  N.Y.,  Macmillan, 
1962. 

6.  Ballistic  Research  Laboratory.  A  Procedure  for  Reading  and  Smooth¬ 
ing  Pressure-Time  Data  from  H.E.  and  Nuclear  Explosions,  by  N.  H. 
Ethridge.  Aberdeen  Proving  Ground,  Md . ,  BRL,  September  1965.  (BRL 
Memorandum  Report  No.  1691.) 

7.  - ,  The  Dependence  of  Blast  on  Ambient  Pressure  and  Temperature, 

by  R.  G.  Sachs.  Aberdeen  Proving  Ground,  Md.,  BRL,  May  1944.  (BRL 
Report  466.) 


19 


UNCLASSIFIED 


UHIUUU 


DOCUMENT  CONTROL  DATA  - R  &  D 


Security  tins  til  ic  at  ion  nl  title,  htnly  nl  otistnn  t  mnl  indeyinf  onimttiti  on  mutt  hr  entered  when  the  uvernll  report  it  <•  Ins  sided) 


i  originating  a  c  t  i  vi  t  v  (Corporate  anther)  p,  report  security  Classification 


Naval  Weapons  Center 
China  Lake,  California  93555 


S  HKPOPT  TITLE 


BLAST  WAVE  YIELDS  FOR  AN  EXPLOSION 


DESCRIPTIVE  NOTES  (Type  nl  report  and  inclusive  dotes! 


Research  Report 


•>  AuTHORiSi  (First  mime ,  middle  initial.  Inst  name) 


Gilbert  F.  Kinney 


UNCLASSIFIED 


October  1970 


8 4t.  CONTRACT  OR  GRANT  NO 


h.  PROJEC  T  NO 


AIR  TASK  A350-350D/216F/0  F17-353-501 


7/1.  TOTAL  NO-  OF  PAGES  I  7b.  NO.  OF  REFS 


|  9/».  ORIGINATOR'S  REPORT  NUMBER(S) 


NWC  TP  5020 


9f\  OTHER  REPORT  NO(SI  (Any  other  numbers  that  may  be  assigned 
this  report) 


10.  DISTRIBUTION  STATEMENT 


This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  the  Naval 
Weapons  Center. 


12.  SPONSORING  MILITARY  ACTIVITY 

Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20360 


13  ABSTRAC  T 


A  blast  wave  overpressure-time  history  contains  voluminous  data  which  can  be 
greatly  simplified  by  converting  them  into  equivalent  yield.  This  calculation 
requires  two  independent  parameters,  usually  chosen  from  measurements  on  peak 
overpressure,  Impulse,  standoff  distance,  or  travel  time.  Representative 
computations  are  described.  Concurrence  among  the  results  of  these  for  a 
particular  explosion  indicates  a  simple  explosion  situation  and  reliable  data; 
divergence  implies  some  effect  such  as  focussed  blast,  or  perhaps  errors  in 
measurement.  All  such  information  is  useful  in  evaluation  of  blast  wave  data, 
and  can  contribute  to  our  understanding  of  explosion  phenomena. 


1473 


(PAGE  I) 


S/N  010  I  •  807-  6  80  t 


UNCLASSIFIED 

Security  Classification 


UNCLASSIFIED 

Security  Classification 


key  wo  r  ds 


Blast  waves 
Explosives 
Explosive  yield 
Instrumentation  for  blast 
Spherical  charges 
Scaling  law 


DD 


FORM 

1  NOV  98 


1473  ( BACK ) 


! 


(PAGE 


2) 


ABSTRACT  CARD 


r 


- r 


o 

nl 

O 

in 

p. 

H 


o 

CN 

o 

m 

Pj 

H 

U 


<13 

c 

8 

o 

c 

qj 


3  .H 
a3  iro 
E  w  qj 


o  a) 
cj  Aj 


wort 

ffl  lAJ  tJ 


a, 

X 

QJ 


QJ  (U 


•H  AJ  G 
AJ  >H 

rt  oo 

S  K 

o  o 

E  *rl 


C 

aj 


l  w 

ai  l  ai  rt 

a  o  u 

3  i-i  43 
qj  a.  ”3  a 
aj  x  c  3 
Aj  ai  rt  w 


C  M-l 
qj  o 


o  aj  c  aj 
C  rt  o  a 

O  r— |  'H  OJ 

u  3  um 
a  rt  iw 
I  3  a) 

•  AJ  AJ 

•a  n  -h  aj 
a)  rt  w  6 
jj  a  o 

T-t  C  W 

w  rt  o 
a  *h  w 
Aj  w  qj 
_  o  o  *H 

«  HJ  H  H 

n.  cl. 
a)  <u  X  6 
Aj  w  <u  -h 
rt  qj 
4=  rt  o 
W  AJ  rH  CJ 

c  a,  c 

o  >AJ  E  rt 

•H  O  i-l  60 
AJ  W  Aj 

rt  w  rt 

AJ  AJ  rt  > 
3  rH  1-i 
CA  3  W  TJ 
B  W  rt 
O  QJ  AJ  •  * 

a  aj  «  rt 

CJ  AJ 
QJ  QJ  l-l  rt 
>  43  rtJ  rtJ 

•H  AJ  c 

AJ  -H  rt 
rt  60  H 
AJ  C  C  43 

c  o  o  rt 
QJ  E  -H  *rl 
W  rt  W  iH 


c 

o 

•H 


E 

QJ 
Aj 
3 

W  AJ 

rt  rt 
aj  3 
R  rH 
rt 

C  > 
•H  rt 

w  C 
Aj  ’H 
O 

Aj  rH 
Aj  3 
QJ  IAJ 
rt 

W  W 

43  W 
Aj  -H 
QJ 

a,  c 
o 

Aj  -H 
O  AJ 
rt 
»  E 

AJ  Aj 
W  O 
rt  4-j 
<H  C 
42  -H 

rtJ  43 
QJ  U 
W  3 
W  W 
3 

O  .H 
O  H 

4J  < 


I 

Aj 

QJ 

TJ 


3 

43  J 
•h  rt 
Aj  C 
aj  aj 
C  E 
o  o 
O  c 

QJ 

C  43 
rt  a. 
u 

c 

•o  o 
S  >H 

rt  w 

o 

*  rH 

rt  CL 
u  X 
rt  «u 
Hd 

4-1 
Q)  O 
> 

rt  oc 

£  c 

•H 

4-»  ’X) 

W  C 
«  to 

iH  4J 
42  W 


rt 

QJ 

rt 

a 

42 

rt 

QJ 

42 

QJ 

QJ 

42 

(0 

3 

Oj 

42 

QJ 

QJ 

(fi 

4-1 

rH 

o 

V4 

•H 

a 

w 

4-1 

rH 

CJ 

rt 

a. 

a 

QJ 

rt 

C 

& 

rt 

o 

4-4 

6 

QJ 

Oj 

C 

Hd 

o 

o 

4-^ 

a 

QJ 

•H 

O 

•H 

OO 

O 

2 

•H 

♦H 

O 

•H 

00 

U 

(0 

M 

U 

•H 

rt 

cn 

4-» 

(0 

u 

rt 

CO 

QJ 

O 

4-1 

o 

rt 

(0 

QJ 

■u 

4J 

rt 

> 

rt 

H 

4-» 

4J 

rt 

> 

rt 

rH 

c 

4-1 

o 

u 

Aj 

rt 

rt 

rH 

rt 

4-1 

4-1 

42 

•H 

QJ 

O 

CJ 

CJ 

4-1 

42 

•H 

QJ 

O 

rt 

> 

O 

QJ 

QJ 

•H 

rt 

> 

HD 

Hd 

43 

rt 

OO 

m 

> 

43 

Hd 

Hd 

Hd 

42 

rt 

OO 

QJ 

o 

S 

c 

AJ 

rt 

QJ 

CJ 

£ 

rt 

QJ 

W 

3 

•H 

Pj 

4-1 

•H 

QJ 

(0 

rt 

•H 

rH 

W 

w 

4-1 

Hd 

H 

rt 

60 

rH 

(0 

w 

■ui 

Hd 

rQ 

rt 

(0 

c 

4-1 

c 

C 

rO 

3 

(0 

c 

rt 

CJ 

rH 

rt 

rt 

U 

c 

O 

o 

rt 

a 

rH 

rt 

rt 

•H 

o 

H 

rH 

4-1 

QJ 

E 

•H 

•H 

0 

rH 

rH 

■U 

rH 

4-1 

< 

42 

(0 

§ 

W 

<3 

W 

rH 

4-1 

< 

42 

(0 

4 


ABSTRACT  CARD 


«^*N 

CO 

P 

3 

1  tH 

3 

•H 

M  >. 

D  rH 

| 

> 

p. 

3  XI 

a  ih 

E 

O 

0 

rH  03 

•H 

1 

o 

• 

• 

•  3  XI 

3  3 

3 

Cp  o 

Q  P  3 

o  in 

P 

00 

W  3  -H 

3 

3 

a 

* 

CTv 

M  O 

CO 

P 

0) 

O  rH 

h  OH 

*H  •* 

3  PS 

H  t— { 

XI  3 

CO 

p 

tH 

p 

in  >>  EH  p 

CO 

03 

XI 

aj 

co  P  B 
<  O  -H 

3 

3 

CJ 

XI 

p 

3 

c 

o 

u 

►J  P  W 

u  in 

x  g  & 

ti 

•H 

rH 

O 

a 

Z  -H  tH  H  3 

0) 

■U 

•H 

o 

O  £  H 

a)  m 

> 

CO 

4J 

•H  03 

o 

rH 

O 

#s 

3  3  3 

tH  P. 

3 

iH 

u 

O  B  3 

> 

p. 

pc 

•H  *H  M 

U  H 

co 

03 

x 

P  P  M 

3  3 

3 

p 

- - 

■*\ 

[d 

03  1 

3  3 

& 

4J 

f 

A 

0  3  3 

rH  Xf 

t 

3 

• 

•H  p  XI 

3  3 

3 

p 

V 

J 

cd 

CM 

rH  3 

>  3 

0 

O 

*H 

Xl  cn  3 

*h  a 

p 

rH 

3  W  3 

3  3 

CO 

(1) 

o 

03 

p.  3  O 

cr  xj 

4J 

CJ 

CM 

u 

P 

3  3 

3 

3 

«  ax: 

•H 

SJ 

3 

p 

CO 

e 

O  P  o 

O 

B 

3 

03 

3 

<N  <D  *H 

H  O 

3 

CO 

M 

H 

J<! 

o  >  x: 

3  S 

P 

*H 

3 

3 

03 

m  o  a 

•H  AJ 

3 

T) 

0) 

*H 

rJ 

CO 

U 

>* 

a  3  3 

E  3 

03 

<H 

03 

H  >  P 

3  3 

3 

IM 

CO 

3 

3 

3  3 

XI  P 

H 

O 

3 

> 

*H 

o  so 

u  -H 

B 

o 

T) 

O  3 
P.  z 

5  §  P  w 

3 

M  O* 

§ 

cd 

^  in  3 

3  3 

p 

4J 

3 

4J 

3  O 

*H  M 

<M 

CO 

Dc 

CO 

• 

*H  c 

4J 

3 

tH 

•  XI  -H 

P  3 

3 

•» 

rH 

rH 

3 

a  B 

3  O 

3 

3 

CO 

CQ 

3 

P.  <!  3 

>  -H 

CO 

CO 

> 

3 

T—\ 

3  P 

o 

-3 

cd 

•H 

o  p 

O  3 

X! 

z 

CN  > 

O  rH 

U 

a 

i  tH 

ia  3ri 

UH 

^  *H  f— i  3  1 

1  •  '  B  t)  8  3 

o  Q  U  0)  o  tn 

£  W  (3  -H  3 

1  cr>  (_|  o  in  w 

rH  Px  U  tH  *H  * 

M  H  £  ffl 

■>  m  w  tH  a  H  i-i 

<U  C/S  P  E  3 

•a  <!  o  in  p 

X  o  ,J  4J  ®  •  a) 

!  P  CJ  (0  X)  E 

1  o  Z  -H  tH  H  3 

o  5  x:  h  <u  >-i 

I  P  *rH  3 

i  «  [!  (II  B  ^  p. 

I  CJ  O  B  3 

l,  2  iH  -H  p  P  P 

i  z  p  v  bo  a  c 

I  CO  I  OJ  111 

*  O  HI  «  H  B) 

|  •  -rl  P  ,0  CO  G 

f  w  h  3  >  a) 

•H  X>  W  C  -H  P- 

i  H  3  II  B  3  m 

B  a  hi  o  o"o 

I  CJ  p  Q)  C 

-ax:  -h 

i  -  o  p  o  o 

I  1)N  dJll  «  p 

1^0  >£  3  5 

I  B  (B  O  S  rl  P 

I  rJ 

i  p*  3  co  B  cn 

B  H  >  U  II  B 

1C  co  co  JZ  P 

••HU  S  x>  P  -H 

l£3  3 

:  U  Z  P  0)  bO  CT 

h-t  cn  3  c  3 

)  co  O  -H  p  ' 


o3 

tH 

•  Xl 

•H 

P 

e 

e 

rH 

rH 

3 

P.  ^ 

B 

3 

o 

3 

3 

03 

PQ 

3 

P.  < 

3 

> 

•H 

m 

V) 

> 

3 

rH 

a 

U 

o 

rH 

cd 

■H 

O 

O 

o 

03 

x: 

3 

55 

CN 

> 

o 

rH 

o 

P- 

t  * 

in  >->  3  >h  I 
C  XI  O  H  B 

•H  rH  CO  »H 

co  -a  co  3 

xi  oj  o  in  « 

C  -H  3  3 
o  w  in  p 

U  -H  *rl  -  3 

h  x  in  in 


tH 

5-4 

cn  tn  P.  H 

P 

in  • 

X2 

0) 

53  ^ 

3 

3  3 

rO 

P 

P  B 

n 

o 

J  p  in  ■ 

Si 

P.  'H 

4J 

u  tn  x 

B 

P  P 

o 

U 

Z  -H  >>rH 

3 

3 

*H 

O 

t)  Z  H  3 

P 

>  rH 

-  C  hi 

o  o  B 


•>0  0) 
••HP 
4-1  H  3 

•h  ^  in 
h  3  in 
to  a  3 

U  P 

r  p, 

•>  O  p 
a)  cn  oj 
yt  o  > 
io  in  o 
►J 

CP  3 
CO  H  > 


3  tH  P* 

3  ^ 

P  p  P  cO 
meet 
3  3  P- 
11  H  O 
XI  3  3  3 
>30 
3  -H  & 
co  3  3  in 

o  cr  -a  p 

3  3  3 


•rl  P  O  3 

£  3  5  p 
S  H  Jj  3 
3 

3  B  V)  3 
P  3  3  3 
3  XS  P  6 
3  P  H 


tH  •  X3 
3  P. 

3  a  < 

3 

•H  O 
t*S  cn 


^'5b 

in  m  o'  o 

3  3  3  P 
O  *H  P  IP 
3  P 

•H  P  3  3 
B  3  O  3 
3  >  iH  in 
rH  3  P  O 
o  o  3  x: 
>  a  h  o 


1  tH 

in  tH  3  rH  I 

3X1  OH  B 

s — \  *H  rH  3  *H 

■  •  •  3  X)  3  3 

i  O  P  p  3  O  V)  » 
I"-  W  3  -H  3  3 

'  <T»  IH  O  4H  V>  P 
H  (u  O  H  -H  -3 


>  p  in  ^  a 
3  cn  p  E 

XI  <  O  -H 


3  3 
in  P 
•H  -  3 
xl  in  w 
H  p  in 

3  3 
P  P 


•  O  U  P  V)  •  3  P. 

p  u  in  T>  B  P 

O  2  -H  >>  rH  3  3 

O  E3  XI  rH  3  p  > 

P  *H  3  O 

«  3  3  3  tH  p. 

CJ  O  B  3  A! 

•  Z  -H  -H  P  P  P  3 

2  P  P  bO  a  a  3 

3  1  3  3  P. 

«  O  3  3  rH  33 
•  -H  p  XI  3  3  3 

4-1  rH  3  >30 

•H  .g  in  3  -H  P. 
i  rH  3  V)  3  3  3  tfl 

3  P.  3  3  0*33  P 

U  p  3  3  3 

•  *rl  3 

-  o  p  o  o  B 

3  CN  3  'H  P  O  3 

^  O  >  Z  3  3  P 

3  CO  O  3  'H  P  3 


3 

a 

0) 

03 

B 

w 

3 

4-1 

H 

> 

U 

3 

3 

3! 

M-4 

a 

03 

cd 

XI 

P 

B 

O 

•H 

y 

3  T3 

p 

•H 

03 

JS 

3 

B 

a 

u 

z 

4J 

co 

bO 

cr 

o 

td 

H-r 

CO 

a 

3 

P 

4J 

03 

o 

•H 

P 

44 

CO 

1  •  rH  3  P 

j  tH  •  X3  -H  P  3  3 

H  3  P.  6  3  0  3 

3  3  P-<S  3  >  -H  « 

3  rH  3  P  O 

•HO  O  O  3  X! 

t«S  CN  >  O  rH  O 


I  to 

<U  I  <U  3 
CJ  O  3 
3  H  X 

i)  at)  o 
1-1  X  3  3 
3  ra  ra  to 
3 

O  3  C  3 

c  ra  o  O 

O  H  -H  t) 
O  3  U  >M 
O  (II  in 
•H  3  0) 

•  3  u 
•a  3  -h  a) 
<u  m  to  6 
o  a  o 
•H  C  to 
3  ra  o 
a  t(  to 
to  3  to  tu 
\  ra  o  o  "H 
ra  >w  h  H 
/  a.  a 
a)  tu  X  E 
3  to  tu  -H 
ra  ra 
X  tu  0) 

to  W  r-t  CJ 

3  ,  x  s 
on  E  til 

•rl  O  *H  60 
3  to  3 
ra  to  tu 
3  3  ra  > 

3  H  -rl 

(X  3  to  *3 
E  to  tu 
O  tu  XI  •- 
o  3  ra  ra 
o  3 
tu  tu  ‘H  ra 
>  x  -a  -a 

'rl  U  3 
XI  "H  tU 

ra  oo  *-< 

3  3  3  43 

3  o  o  ra 

tU  E  'rl  *H 

to  ra  to  rH 


•  i 

XI  1-1 

C  txi  tu 
tu  O  T3 
S  3 
tu  3  3 

£  ° 

3  -H  3 
to  xi  3 
ra  ra  o 
<u  3 
B  H  O 

ra  xi 

3  3* 

•H  <U  (U 
XI 

to  3  3 

It  rl 
O  -H 
3  H  3 
3  3  Xi 

tin  3 
tu  o 
to  to  a 
(X  3 
ra  3 
X  to  ra 
3  3  u 
tu 

a  b  ra 

O  3 
3  ’rl  t0 

O  3 

ra  r 
r  6  ra 

3  C  3 
to  o  ra 
ra  «3  tj 

rH  C 

X  ’H  tu 
> 

•a  x  ra 
tu  u  s 
to  3 
to  to  3 


I  to 

tu  i  tu  ra 

0  0  3 
3  rl  43 

tu  tx  *a  o 
3X33 
3  tu  ra  to 
3 

0  3  3  3 

3  ra  o  o 

OH'H  (U 
(J  3  3  13 

O  ra  txt 
-H  3  tu 

♦  3  3 

Tl  3  rl  U 
tu  ra  to  £ 


3  ra  o 
O  ’H  to 
to  3  to  tu 
(U  O  O  ’H 

V3  C3  r-l  r-l 

a  a 
7tU  tU  X  E 
3  to  tU  ’H 
ra  tu 
X  <u  <u 

to  3  r-l  O 
3  as 
o  n  E  j 


3  3  ra  > 

3  rH  *rl 

a  3  to  ra 
E  to  ra 

o  a)  3  .r 

o  3  ra  ra 

o  3 

tu  (U  rl  ra 

>  X  XI  X 
•rl  3  a 
3  ’rl  (U 
ra  60  r-l 

3  3  3  43 
3  o  o  ra 

tU  E  ’rl  ’r-l 
to  ra  to  r-l 


I  to 
tu  i  tu  ra 
0  0  3 
3  r-l  X 
tu  fx  ra  o 
3  X  3  3 
3  tu  ra  to 
3 

0  3  3  3 

3  ra  o  o 
O  H  rl  IU 
U  3  3  >3 

o  ra  tx 
•h  3  tu 

•  3  4J 

ra  3  ’H  01 
tu  ra  to  E 
43  a  o 
■h  s  to 
3  ra  o 

O  rH  to 
to  3  to  (U 
\  tu  O  O  ’rl 


tu  tu  &  E 
3  to  tu  ’rl 
ra  tu 

X  tu  (U  , 


O  n  e  tu 

•rl  O  ’rl  60 


3  3  ra  > 

3  r-l  r-l 

X  3  to  X) 
E  to  tu 
o  in  u  •» 
o  3  ra  ra 

O  3 

ra  tu  -h  ra 
?  3  ra  ra 

•H  3  C 
3  -rl  tU 

ttl  60  r-l 

3  3  3  43 

s  o  o  ra 

(U  6  rl  'rl 
CO  ra  to  r-l 


